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A general theory for the condensation of strongly bound quartets in infinite nuclear 
matter is presented. Critical temperatures for symmetric and asymmetric nuclear matter 
are evaluated. A fully nonlinear theory for the quartet order parameter, based on an analogy 
of the Gorkov approach to pairing, is presented and solved. The strong qualitative difference 
^ ■ with pairing is pointed out. 



O ■ §1. Introduction 

(N 

Since about one decade, there has been a renewed and strong interest in gas 
like a particle states in nuclei, as, e.g. the famous Hoyle state in 12 C. This stems 
essentially from the fact that in 2001 the prediction was made that besides the 
Hoyle state with three a's, other similar states with more a's can exist and that 
those states can be considered as being of a particle condensate typeP Of course, 
in finite systems, like nuclei, neither pairing nor quartetting can show phase locking 
as in macroscopic systems. They can be considered as precursors to condensation 
phenomena in macroscopic systems, like neutron or nuclear matter. The proposed 
THSR wave function 1 ' for the description of the a gas states is very successful for 
nuclei. It is a number conserving ansatz, in spirit similar to a number projected BCS 
wave function, which cannot be extended to the infinite matter case. We, therefore, 
\Q \ will discuss in this article how a particle condensation can be described for the case of 

infinite nuclear matter. This may be relevant for compact stellar objects. However, 
the general theory may also serve for other systems where quartet correlations are 
important as, e.g., in a gas of excitons where bi-excitons may be formed. Also in the 
physics of cold atoms the possibility of trapping four different fermions and, thus, 
of quartet condensation, may exist in the future. The paper is organised as follows. 
In £}2j we will describe how to determine the critical temperature in symmetric and 
asymmetric matter. In [j3j a general formula for the quartet order parameter will be 
presented and applied to symmetric nuclear matter using a Gorkov type of formalism. 
Finally in SJHwe give some perspectives and conclude. 

§2. Critical temperature of quartetting and a particle condensation 

In this section, we will discuss the critical temperature of quartet condensation in 
symmetric and asymmetric nuclear matter. A comparison with the pairing case will 
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also be given. For this, we start out to establish an in medium four body equation. 

In medium four particle correlation for example appears if one adds an a-particle 
on top of a double magic nucleus such as 208 Pb, or in semiconductors where in the 
gas of excitons, i.e. p — h bound states, may appear bi-excitons, i.e. bound state of 
two p — h pairs. The effective wave equation contains, in mean field approximation, 
the Hartree-Fock self-energy shift of the single-particle energies as well as the Pauli 
blocking of the interaction. We consider for the derivation the Equation of Motion 
Method (EMM), well known from the Random Phase Approximation (RPA)PJ 

([SQ,[H,SQl t }]) = E p ([SQ,Qi]), (2-1) 

where 

G« = V II E ^(1234)4444, (2-2) 

' 1234 

is the quartet creator. Of course, we also could consider analogously a quartet 
destructor. Considering a fermion Hamiltonian with two body interaction, we eval- 
uate the double commutator in ()2Tp and perform the expectation value with the 
Hartree-Fock ground state, which is, as we know, the standard procedure to obtain 
RPA equations in linearised form. The result iJS® 

{E p - e x - e 2 - e 3 - e 4 )^4, P (1234) 
= [1 - h ~ /2]^i2i'2'V'4,p(l / 2 / 34) + permutations, (2-3) 

where ei is the single particle energy, i.e. the kinetic energy plus HF shift. A similar 
equation can be obtained for three particles, i.e. A = 3. This equation has exactly the 
same structure as a free four body Schrodinger equation with the only difference that 
the two body matrix elements Viji>j> are premultiplied with the phase space factor 
fifj — fifj = 1 — fi — fj where the /, are the Fermi-Dirac functions, and /i = 1 — /j. 
The Fermi function is equal to the step function at T = 0, i.e. fi = 0(n — e^), but 
immediately generalisable to finite temperature with /j = [1 + e^ ei ~ fl ^ T ]~ 1 , a result 
which could have been derived employing Matsubara Green functions. In (|2-3p the 
€i are as usual the HF single particle energies and t^yy = (ij\V\i'f) — (ij\V\j'i') 
is the antisymmetrized matrix element of the two body interaction. We realise that 
the phase space factor in (|2-3p is exactly the same as in the two particle RPAP It 
turns out to be a general rule in deriving higher pp or ph RPA's that the pp and ph 
matrix elements are always premultiplied with the same phase space factors as in 
the corresponding standard pp or ph RPA's. 

We are interested in an example of nuclear physics where the a-particle con- 
stitutes a particularly strongly bound cluster of four nucleons. One can ask the 
question how, for a fixed temperature, the binding energy of the a-particle varies 
with increasing temperature. 

The effective wave equation has been solved using separable potentials for A = 2 
by integration. For A = 3, 4 we can use a Faddeev approach^ The shifts of binding 
energy can also be calculated approximately via perturbation theory. In Fig. Q] 
we show the shift of the binding energy of the light clusters (deuteron (d), triton 
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(t), helion (h) and a) in symmetric nuclear matter as a function of density for 
temperature T = 10 MeVp* 

It is found that the cluster bind- 
ing energy decreases with increasing 
density. Finally, at the Mott den- 
sity PA°n*p( T ) the bound state is dis- 
solved. The clusters are not present at 
higher densities, merging into the nu- 
cleonic medium. It is found that the 
a particle already dissolves at a den- 
sity p^ ott fa po/W, see Fig. Q] For a 
given cluster type characterized by A, n, 
we can also introduce the Mott momen- 
tum P^ I ° tt (/9, T) in terms of the ambient 
temperature T and nucleon density p, 
such that the bound states exist only for 
P > Pf° u (p,T). We do not present an 
example here, but it is intuitively clear 
that a cluster with high c.o.m. momen- 
tum with respect to the medium is less 
affected by the Pauli principle than a 
cluster at rest. 

Since Bose condensation only is of relevance for deuterons (d) and a's, and the 
fraction of d, tritons (t) and helions (h) becomes low compared with that of a's with 
increasing density, we can neglect the contribution of the latters to the equation 
of state. Consequently, if we further neglect the contribution of the four-particle 
scattering phase shifts in the different channels, we can now construct an equation 



Fig. 1 



Shift of binding energy of the 
light clusters (d - dash dotted, t/h - dot- 
ted, and a - dashed: perturbation the- 
ory, full line: non-perturbative Faddeev- 
Yakubovski equation) in symmetric nu- 
clear matter as a function of density for 
given temperature T = 10 MeV. 



of state p(T, p) = p tTee (T,p) + p 



.bound, d 



(T,p)+P 



bound, a 



(T, p) such that a-particles 



determine the behavior of symmetric nuclear matter at densities below p^ ott and 
temperatures below the binding energy per nucleon of the a-particle. The formation 
of deuteron clusters alone gives an incorrect description because the deuteron binding 
energy is small, and, thus, the abundance of d-clusters is small compared with that 
of a-clusters. In the low density region of the phase diagram, a-matter emerges as 
an adequate model for describing the nuclear-matter equation of state. 

With increasing density, the medium modifications - especially Pauli blocking - 
will lead to a deviation of the critical temperature T c (p) from that of an ideal Bose 
gas of a-particles (the analogous situation holds for deuteron clusters, i.e., in the 
isospin-singlet channel) . 

Symmetric nuclear matter is characterized by the equality of the proton and 
neutron chemical potentials, i.e., p p = p n = p. Then an extended Thouless condition 
based on the relation for the four-body T-matrix (in principle equivalent to (|2-3p at 
eigenvalue Ap) 



T 4 (1234, 1"2"3"4", Ap) = ^ 



1'2'3'4' 



4/x - ei - e 2 - e 3 - e 4 



<5(3,3')5(4,4 / ) +cycl. 
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xT 4 (l , 2'3'4 / , 1"2"3"4", 4/x) (2-4) 

serves to determine the onset of Bose condensation of a- like clusters, noting that 
the existence of a solution of this relation signals a divergence of the four-particle 
correlation function. An approximate solution has been obtained by a variational 
approach, in which the wave function is taken as Gaussian incorporating the correct 
solution for the two-particle problem.® 
On the other hand, Eq. (|2-4|) . re- 
spectively (|2-3p at eigenvalue 4//, has 
also been solved numerically exactly by 
the Faddeev-Yakubovsky method em- 
ploying the Mamiet-Tjon force.®' The 
results for the critical temperature of a- 
condensation is presented in Fig. [2] as 
a function of the chemical potential //. 
The exact solution could only be ob- 
tained for negative /i, i.e. when there 
exists a bound cluster. It is, therefore, 
important to try yet another approxi- 
mate solution of the in-medium four- 
body equation. Since the a-particle is 
strongly bound, we make a momentum 
projected mean field ansatz for the quar- 
tet wave function^ « 

4-1234 = (27r) 3 ^ 3 )(k! + k 2 + k 3 + k 4 ) 
4 

x]J l p(k i )x ST , (2-5) 

i=\ 

where \ ST is the spin-isospin function 
which we suppose to be the one of a 
scalar (S = T = 0). We will not fur- 
ther mention it from now on. We work 
in momentum space and <£>(k) is the as-yet unknown single particle OS wave function. 
In position space, this leads to the usual formulsP ^1234 — > J d 3 RY\f =1 </5(rj — R) 
where </3(rj) is the Fourier transform of y?(kj). If we take for ty?(kj) a Gaussian 
shape, this gives: ^1234 — > exp[— c ^] 1<i<fc<4 (rj — r^) 2 ] which is the translationally 
invariant ansatz often used to describe a-clusters in nuclei. For instance, it is also 
employed in the a-particle condensate wave function of Tohsaki, Horiuchi, Schuck, 
Ropke (THSR) in Ref.E]). 

Inserting the ansatz ()2-5p into (|2-3p and integrating over superfluous variables, 
or minimizing the energy, we arrive at a Hartree-Fock type of equation for the single 
particle OS wave function ip(k) = </?(|k|) which can be solved. However, for a general 
two body force ^ klka k / k ^, the equation to be solved is still rather complicated. We, 
therefore, proceed to the last simplification and replace the two body force by a 
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Fig. 2. Critical temperature of alpha and 
deuteron condensations as functions of 
(a) chemical potential and (b) density 
of free nucleon Crosses (x) corre- 

spond to calculation of Eq. (|2-3p with the 
Malfliet-Tjon interaction (MT I-III) using 
the Faddeev-Yakubovski method. 
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unique separable one, that is 

*Wm = Xe~ k2 / k oe- k ' 2 / k h2nf5^(K - K% (2-6) 

where k = (k x - k 2 )/2, k' = (ki - k' 2 )/2, K = k x + k 2 , and K' = k[ + k 2 . This 
means that we take a spin-isospin averaged two body interaction and disregard that 
in principle the force may be somewhat different in the S, T = 0, 1 or 1,0 channels. 
It is important to remark that for a mean field solution the interaction only can be 
an effective one, very different from a bare nucleon-nucleon force. This is contrary 
to the usual gap equation for pairs, to be considered below, where, at least in the 
nuclear context, a bare force can be used as a reasonable first approximation. 

We are now ready to study the solution of (|2-3p and (|2-4p for the critical temper- 
ature T", defined by the point where the eigenvalue equals 4/i. For later comparison, 
the deuteron (pair) wave function at the critical temperature is also deduced from 
the Thouless criterion for the onset of pairing with (|2-6p to be 

Hk) = _ 1 - 2 /( £ ( fc )) Ae -^ / ^e- k2 l k h(k% (2-7) 



k 2 /m-2^ J (2tt) 

where 4>(k) is the relative wave function of two particles given by $12 — > 0(| kl 7 2 |) 
«J( 3 ) (ki + k 2 ), ande(£;) = k 2 /{2 m)—fj,. We also neglected the momentum dependence 
of the Hartree-Fock mean field shift in Eq. (|2-7p . It, therefore, can be incorporated 
into the chemical potential \i. With Eq. (|2-7p . the critical temperature of pair con- 
densation is obtained from the following equation: 

f d?k l-2/(e(fe)) 2fc a /fc g 

1 -~ X J (2tt)3 k y m -2v e ■ (2 8) 

In order to determine the critical temperature for a-particle condensation, we 
have to adjust the temperature so that the eigenvalue of (|2-3p and (|2-4p equals 4/x. 
The result is shown in Fig. [21(a) . In order to get an idea how this converts into a 
density dependence, we use for the moment the free gas relation between the density 
riA°) of uncorrelated nucleons and the chemical potential 

n(0) = 4 / (0 /(£) - (2 ' 9) 

We are well aware of the fact that this is a relatively gross simplification, for instance 
at the lowest densities, and we intend to generalize our theory in the future so that 
correlations are included into the density. This may be done along the work of 
Nozieres and Schmitt-RinkA ! ' The two open constants A and ko in Eq. (|2-6p are 
determined so that binding energy (—28.3 MeV) and radius (1.71 fm) of the free 
{ft = 0) a-particle come out right. The adjusted parameter values are: A = —992 
MeV fm 3 , and ko = 1.43 fm -1 . The results of the calculation are shown in Fig. [2J 

In Fig. [21 the maximum of critical temperature T c a max is at \i = 5.5 MeV, and 
the a-condensation can exist up to /x max = 11 MeV. It is very remarkable that the 
results obtained with (|2-5|) for T" very well agree with the exact solution of (|2-3p and 
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(|2-4p using the Malfliet-Tjon interaction (MT I-III) 61 with the Faddeev-Yakubovski 
method also shown by crosses in Fig. [2] (the numerical solution only could be obtained 
for negative values of fi). This indicates that T" is essentially determined by the 
Pauli blocking factors. 

In Fig. [2] we also show the critical temperature for deuteron condensation derived 
from Eq. (j2-8|) . In this case, the bare force is adjusted with A = —1305 MeV fm 3 
and fco = 1-46 fm -1 to get experimental energy (—2.2 MeV) and radius (1.95 fm) 
of the deuteron. It is seen that at higher densities deuteron condensation wins over 
the one of a-particles. The latter breaks down rather abruptly at a critical positive 
value of the chemical potential. Roughly speaking, this corresponds to the point 
where the a-particles start to overlap substantially. This behavior stems from the 
fact that Fermi-Dirac distributions in the four body case, see (|2-4p . can never become 
step-like, as in the two body case, even not at zero temperature, since the pairs in 
an a-particle are always in motion. Therefore, no threshold effect occurs as with 
pairing for Cooper pairs at rest. As a consequence, a-condensation generally only 
exists as a BEC phase and the weak coupling regime is absent. 

An important consequence of this study is that at the lowest temperatures, 
Bose-Einstein condensation occurs for a particles rather than for deuterons. As the 
density increases within the low-temperature regime, the chemical potential [i first 
reaches —7 MeV, where the a's Bose-condense. By contrast, Bose condensation of 
deuterons would not occur until \x rises to —1.1 MeV. 

The "quartetting" transition temperature sharply drops as the rising density 
approaches the critical Mott value at which the four-body bound states disappear. 
At that point, pair formation in the isospin-singlet deuteron-like channel comes into 
play, and a deuteron condensate will exist below the critical temperature for BCS 
pairing up to densities above the nuclear-matter saturation density po> as described 
in the previous Section. Of course, also the well known and studied isovector n-n and 
p-p pairing develops. Usually, p-n pairing in the isoscalar channel is stronger. How- 
ever, the competition between the different pairing channels is not well known. The 
critical density at which the a condensate disappears is estimated to be po/3. There- 
fore, a-particle condensation primarily only exists in the Bose-Einstein-Condensed 
(BEC) phase and there does not seem to exist a phase where the quartets acquire 
a large extension as Cooper pairs do in the weak coupling regime. However, the 
variational approaches of Ref. [3]) and of Eq. (|2-5p on which this conclusion is based 
represent only first attempts at the description of the transition from quartetting 
to pairing. The detailed nature of this fascinating transition remains to be clari- 
fied. Many different questions arise in relation to the possible physical occurrence 
and experimental manifestations of quartetting: Can we observe the hypothetical "a 
condensate" in nature? What about thermodynamic stability? What happens with 
quartetting in asymmetric nuclear matter? Are more complex quantum condensates 
possible? What is their relevance for finite nuclei? As discussed, the special type 
of microscopic quantum correlations associated with quartetting may be important 
in nuclei, its role in these finite inhomogeneous systems being similar to that of 
pairingpJ 1 

On the other hand, if at all, a-condensation in compact star occurs at strongly 
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asymmetric matter. It is, therefore, important to generalize the above study for 
symmetric nuclear matter to the asymmetric case. This can be done straightfor- 
wardly, again using our momentum projected mean field ansatz (|2-5|) generalized to 
the asymmetric case. This implies to introduce two chemical potentials, one for neu- 
trons and for protons. We also have to distinguish two single particle wave functions 
in our product ansatz which now reads® 

■01234 -> Vp{ki)(p P (k2)(Pn(k 3 )iPn(k / t)xo(2'K) 3 5(k 1 + k 2 + k 3 + k 4 ) (2-10) 

where ty T {ki) = yv(|kj|) is the s-wave 
single particle wave functions for pro- 
tons (r = p) and neutrons (r = n), re- 
spectively, xo is the spin-isospin singlet 
wave function. This now leads to two 
coupled equations of the Hartree-Fock 
type for (p n and ip p . For the force we 
use the same as in the symmetric case. 

Fig. E^a) shows the critical temper- 
ature of a condensation as a function 
of the total chemical potential /U to tai = 
AV + fJ"n- We see that T c decreases as the 
asymmetry, given by the parameter 
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increases. This is in analogy with the 
deuteron case (also shown) which al- 
ready had been treated in Refs. fTU |) .[TT ]) . 
On the other hand, in Fig. [3^b), it is 
also interesting to show T c as a function 
of the free density which is 
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Fig. 3. Critical temperature as functions of 
the total chemical potential /itotai — + 
fi n (top) and the total free density n tota i 
(bottom) Thick (thin) lines are for a- 
particle (deuteron). Solid, dashed, and 
dotted lines are respectively for 5 = 0.0, 
S = 0.5, and S — 0.9, where the density 
ratio S is in Eq. (|2- 1 1 [I . 



where the factor two in front of the in- 
tegral comes from the spin degeneracy, 
and fp, n (k) = [l + exp(A; 2 /2m-/ip iri )] _1 . 
It should be emphasized, however, that 

in the above relation between density and chemical potential, the free gas relation 
is used and correlations in the density have been neglected. In this sense the de- 
pendence of T c on density only is indicative, more valid at the higher density side. 
The very low density part where the correlations play a more important role, will be 
treated in the future. For instance, it will be important to recover the the critical 
temperature corresponding to an ideal Bose gas in that limit. Techniques similar to 
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the one of Nozieres and Schmitt-Rink may be employed.^ It should, however, be 
stressed that the dependence of T c on the chemical potential as in Fig. [3]^ a), stays 
unaltered. It is only the relation between the chemical potential and the (correlated) 
density which changes. 

The fact that for more asymmetric matter the transition temperature decreases, 
is natural, since as the Fermi levels become more and more unequal, the proton- 
neutron correlations will be suppressed. For small <5's, i.e., close to the symmetric 
case, a condensation (quartetting) breaks down at smaller density (smaller chemical 
potential) than deuteron condensation (pairing). This effect has already been dis- 
cussed above for symmetric nuclear matter ®. For <5's close to one, i.e. strong asym- 
metries, the behavior is opposite, i.e., deuteron condensation breaks down at smaller 
densities than a condensation, because the small binding energy of the deuteron can 
not compensate the difference of the chemical potentials. 

More precisely, for small <5's, the deuteron with zero center of mass momentum 
is only weakly influenced by the density or the total chemical potential as can be 
seen in Fig. [3l However, as 5 increases, the different chemical potentials for protons 
and neutrons very much hinders the formation of proton-neutron Cooper pairs in the 
isoscalar channel for rather obvious reasons. The point to make here is that because 
of the much stronger binding per particle of the a-particle, the latter is much less 
influenced by the increasing difference of the chemical potentials. For the strong 
asymmetry 5 = 0.9 in Fig. [3] then finally a-particle condensation can exist up to 
^totai = 0.02 fm -3 (/x tota i = 9.3 MeV), while the deuteron condensation exists only 
up to ntotai = 0.005 fm" 3 (/x total = 6.0 MeV). 

Overall, the behavior of T c is more or less as expected. We should, however, 
remark that the critical temperature for a-particle condensation stays quite high, 
even for the strongest asymmetry considered here, namely 5 = 0.9. This may be 
of importance for the possibility of a-particle condensation in neutron stars and 
super novae explosions. 

In conclusion the a-particle (quartet) condensation was investigated in homoge- 
neous symmetric nuclear matter as well as in asymmetric nuclear matter. We found 
that the critical density at which the a-particle condensate appears is estimated to 
be around po/3 in the symmetric nuclear matter, and the a-particle condensation 
can occur only at low density. This result is consistent with the fact that the Hoyle 
state (0^~) of 12 C also has a very low density p ~ Po/3- On the other hand, in 
asymmetric nuclear matter, the critical temperature T c for the a-particle condensa- 
tion was found to decrease with increasing asymmetry. However, T c stays relatively 
high for very strong asymmetries, a fact of importance in the astrophysical context. 
The asymmetry affects deuteron pairing more strongly than a-particle condensation. 
Therefore, at high asymmetries, if at all, a-particle condensate seems to dominate 
over pairing at all possible densities. 

§3. 'Gap' equation for the quartet order parameter 

For macroscopic a condensation it is, of course, not conceivable to work with a 
number projected a particle condensate wave function as we did when in finite nuclei 
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only a couple of a particles were present P We rather have to develop an analogous 
procedure to BCS theory but generalized for quartets. In principle a number non- 
conserving wave function of the type \a) = exp[^ 1234 z\234cf c£ c% ]|vac) would be 
the ideal generalization of the BCS wave function for the case of quartets. However, 
unfortunately, it is unknown so far (see, however, Ref. I12p ) how to treat such a 
complicated many body wave function mathematically in a reasonable way. So, we 
rather attack the problem from the other end, that is with a Gorkov type of approach, 
well known from pairing but here extended to the quartet case. Since, naturally, the 
formalism is complicated, we only will outline the main ideas and refer for details to 
the literature. 

Actually one part of the problem is written down easily. Let us guide from a 
particular form of the gap equation in the case of pairing. We have at zero temper- 
ature 

(ei + e 2 )Kl2 + (1 - Pi ~ P2)^ ^2 Vl21'2'«l'2' = 2/XK12, (3d) 

1'2' 



where ku = {c\C2) is the pairing tensor, p, = (cfci) = \ 



1 



with Ei 



i - p) 2 + Af the quasi-particle energy are the BCS occupation numbers, and 

^i2i'2' denotes the antisymmetrized matrix element of the two-body interaction. The 
Si are the usual mean field energies. Equation (|3-ip is equivalent to the usual gap 
equation in the case of zero total momentum and opposite spin, i.e. in short hand: 
2 = 1 where the bar stands for 'time reversed conjugate'. The extension of (|3- 1[) to 
the quartet case is formally written down without problem and can be derived with 
EMM, slightly extending (|2-3p to correlated occupation numbers at T=0. 



(£1234 - 4^)^1234 = (1 ~ Pi ~ P?) ^ X] ^121'2'«1'2'34 

1'2' 

+ (1 - pi - p 3 )- ^ Vi3i'3'«l'23'4 + a 11 permutations. (3-2) 

1'3' 

with K1234 = (C1C2C3C4) the quartet order parameter. This is formally the same equa- 
tion as in Eq. (|2-3p with, however, the Fermi-Dirac occupation numbers replaced by 
the zero temperature quartet correlated single particle occupation numbers, sim- 
ilar to the BCS case. For the quartet case, the crux lies in the determination 
of those occupation numbers. Let us again be guided by BCS theory or rather 
by the equivalent Gorkov approachPS In the latter, there are two coupled equa- 
tions, one for the normal single particle Green's function (GF) and the other for 
the anomalous GF. Eliminating the one for the anomalous GF in inserting it into 
the first equation leads to a Dyson equation with a single particle mass operator,^ 




Fig. 4. Graphic representation of the BCS 
mass operator in Eq. (|3-3p 
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12 



34K43 



(3-4) 



31 




This can be graphically represented as 

in Fig. H] where (cc) stands for the order parameter k± 2 and the dot for the two body 
interaction. 

The generalization to the quartet case is considerably more complicated but 
schematically the corresponding mass operator in the single particle Dyson equation 
can be represented graphically as in Fig. [5] with the quartet order parameter (cccc). 
This mass operator can again be de- 
rived employing the EMM. However, the 
details are rather lengthy and involved 
and, thus, we refrain from presenting 
this here, see.® Put aside this diffi- 
culty to derive a manageable expression 
for this 'quartet' s.p. mass operator, 
what immediately strikes is that instead 
of only one 'backward going line' with 
(— p, —a) as in the pairing case, we now 
have three backwards going lines. As a 

consequence, the three momenta ki, k2, k3 in these lines are only constrained so 
that their sum be equal to ki + k 2 + k 3 = — p in order that the total momentum 
of the order parameter be zero and, thus, the remaining freedom has to be summed 
over. This is in strong contrast to the pairing case where the single backward go- 
ing line is constrained by momentum conservation to have momentum —p. So, no 
internal summation occurs in the mass operator belonging to pairing. The conse- 
quence of this additional momentum summation in the mass operator for quartetting 
leads with respect to pairing to a completely different analytic structure of the mass 
operator in case of quartetting. This is best studied with the so-called three hole 
level density g 3 h{u) which is related to the imaginary part of the three hole Green's 
function G 3h (k 1 ,k 2 ,k 3 ;uj) = (/1/2/3 + hhh)/{u + £123) with Ei 23 = £1 + e 2 + £3 
figuring in the mass operator, see Fig. [5] (we remind: / = ! — /) 



Fig. 5. Graphical representation of the ap- 
proximate a-BEC mass operator M quartct 
of Eq. (323). 



d 3 k t d 3 k 2 d 3 k 3 



(2tt) 3 (2vr) 3 (2vr) 3 
d 3 h d 3 k 2 d 3 k 3 



ImG^ 3h \k u k 2 ,k 3 ;uj + i7 ] ) 



(2tt) 3 (2tt) 3 (2tt) 3 



(/1/2/3 + fif2h)n8(u> + E X + e 2 + e 3 ). (3-5) 



In Fig. [6] we show the level density at zero temperature (/(w) = 9(—uj)), where 
it is calculated with the proton mass m = 938.27 MeVP Two cases have to be 
considered, chemical potential fi positive or negative. In the latter case we have 
binding of the quartet. Let us first discuss the case \x > 0. We remark that in this 
case, the 3h level density goes through zero at uj = 0, i.e., since we are measuring 
energies with respect to the chemical potential fi, just in the region where the quartet 



Quartetting in Nuclear Matter 



11 



correlations should appear. This is at strong variance with the pairing case where 
the lh level density, g lh (u) = J $fa(f k + f k )5(u> + e(k)) = J ^S(u; + e(k)), 
does not feel any influence from the medium and, therefore, the corresponding level 
density varies (neglecting the mean field for the sake of the argument) like in free 
space with the square root of energy. In particular, this means that the level density 
is finite at the Fermi level. This is a dramatic difference with the quartet case 
and explains why Cooper pairs can strongly overlap whereas for quartets this is 
impossible as we will see below. We also would like to point out that the 3h level 
density is just the mirror to the 3p level density which has been discussed in Ref.ll4p. 



For the case where [i < where 
anyway the fi's are zero at T=0, there 
is nothing very special, besides the fact 
that the three hole level density only is 
non-vanishing for negative values of ui 
and that the upper boundary is given 
by uj = 3/x. Therefore, the level density 
of Eq. (|3-5p is zero for u) > 3/i. There- 
fore, in the BEC regime (fj, < 0), there is 
no marked difference between the pair- 
ing and quartetting cases. 

With these preliminary but crucial 
considerations we now pass to the eval- 
uation of the s.p. mass operator with 
quartet condensation. Its expression 
corresponding to Fig. [6] can be shown 
to be of the following form 



234 



^1234(/2/ 3 /4 + f2hh)A 



1234 



W + ^234 



(3-6) 



with 



^1234 = ^ ^^12,1 / 2'^33'^44'^1'2'3'4 / 

(3-7) 



1'2'3'4' 



Again, comparing the quartet s.p. mass 
operator (|3-6p with the pairing one (|3-3p , 
we notice the presence of the phase 
space factors in the former case while 
in Eq. (|3-3p they are absent. As already 
indicated above, this fact implies in the 




Fit 
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. 6. 3h level densities denned in Eq. (|3-5[) 
for various values of the chemical potential 



fj, at a zero temperature 
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quartet case that only the Bose-Einstein 
condensation phase is born out whereas 
a 'BCS phase' (long coherence length) is 

absent, since the three hole level density is zero at 3fi due to the non vanishing phase 
space factors. 

The complexity of the calculation in Eq. (|3-6|) is much reduced using for the 
order parameter (cccc) our mean field ansatz projected on zero total momentum, as 
it was already very successfully employed with Eq. (|2-5p . 

(C1C2C3C4} ->■ 0k 1 k 2 ,k 3 k 4 Xo, 

0k lk2 ,k3k4 = ^(|k 1 |) V 9(|k 2 |)^(|k 3 |)^(|k4|)(27r) 3 5(k 1 + k 2 + k 3 + k 4 ), (3-8) 

where xo is the spin-isospin singlet wave function. It should be pointed out that 
this product ansatz with four identical OS" single particle wave functions is typical 
for a ground state configuration of the a particle. Excited configurations with wave 
functions of higher nodal structures may eventually be envisaged for other physical 
situations. We also would like to mention that the momentum conserving 5 function 
induces strong correlations among the four particles and (|3-8p is, therefore, a rather 
non trivial variational wave function. 

For the two-body interaction of Vk 1 k2,k 3 k4 m Eq. (|3-7p . we employ the same 
separable form (|2-6p as done already for the quartet critical temperature. 




-250 -200 -150 -100 -50 -250 -200 -150 -100 -50 



w (MeV) uj (MeV) 

Fig. 7. -ImA/ quartct (fci,aj + irj) in Eq. (EHJ as a function of u for p = -4.9MeV (left) and for 
fi = 0.55MeV (right) at zero temperature!^* 

At first let us mention that in this pilot application of our selfconsistent quartet 
theory, we only will consider the zero temperature case. As a definite physical 
example, we will treat the case of nuclear physics with the particularly strongly 
bound quartet, the a particle. It should be pointed out, however, that if scaled 
appropriately all energies and lengths can be transformed to other physical systems. 
For the nuclear case it is convenient to measure energies in Fermi energies £f = 
35 MeV and lengths in inverse Fermi momentum k]? 1 = 1.35 -1 fm. 

We are now in a position to solve, as in the BCS case, the two coupled equations 
(|3Tp for the quartet order parameter and the single particle occupation numbers 
from the single particle Dyson equation with single particle self energy (|3-6p with 
(|3-8p selfconsistently. The single particle wave functions and occupation numbers 
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k (fm- 1 ) r (fm) k (fm- 1 ) 



Fig. 8. Single particle wave function (p(k) in fc-space (left), for r-space <p(r) (middle), and occupa- 
tion numbers (right) at /i = —5.26 (top), —1.63 (middle) and 0.55 (bottom). The r-space wave 
function <fi(r) is derived from the Fourier transform of tp(k) by <p(r) = f d 3 fce lkr ip(k) / (2tt) s . The 
dashed line in the left panels correspond to the Gaussian with same norm and r.m.s. momentum 
as <p(fc)P 



obtained from the above cycle are shown in Fig. [8l We also insert the Gaussian wave 
function with same r.m.s. momentum as the single particle wave function in the left 
figures in Fig. [8) As shown in Fig. [8l the single particle wave function is sharper 
than a Gaussian. 

We could not obtain a convergent solution for [i > 0.55 MeV. This difficulty has 
precisely it's origin in the fact that the three hole level density goes through zero at 
3/U > 0, just where the four body correlations should build up, as this was discussed 
above. In the r.h.s. panels of Fig. [8] we also show the corresponding occupation 
numbers. We see that they are very small. However, they increase for increasing 
values of the chemical potential. For [i = 0.55 MeV the maximum of the occupation 
still only attains 0.35 what is far away from the saturation value of one. What really 
happens for larger values of the chemical potential, remains unclear. 

The situation in the quartet case is also in so far much different, as the 3/i 
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Fig. 9. — 2ImG(A;, u>+irf) in Eq. (|3-6|l as function of ui for /i 
(bottom) at a zero temperature.^ 
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Green's function produces a considerable imaginary part of the mass operator. 

Figure [7] shows the imaginary part of the approximate quartet mass operator 
of Eq. (|3-6p for fi < and /i > 0. These large values of the damping rate imply a 
strong violation of the quasiparticle picture. In Fig. Owe show the spectral function 
of the single particle GF. Contrary to the pairing case with its sharp quasiparticle 
pole, we here only find a very broad distribution, implying that the quasiparticle 
picture is completely destroyed. How to formulate a theory which goes continuously 
from the quartet case into the pairing case is, as mentioned, an open question. One 
solution could be to start right from the beginning with an in medium four body 
equation which contains a supernuid phase. When the quartet phase disappears, the 
supernuid phase may remain. Such investigations shall be done in the future. 



§4. Discussion, perspective, and conculision 



The considerations in £J3] ar e adequate for a situation where the uncorrelated 
Fermi gas directly goes over into the a-particle condensed phase. However, in reality, 
in the low density phase of nuclear matter, the process may go via bound states of 
tritons or helions. One may imagine a mixture of deuterons, tritons, helions, and 
nucleons. The capture of a further nucleon of the triples can lead to a condensation. 
The mass operator of such a process is depicted in Fig. [TUTa). Since the triples now 
have a definite c.o.m. momentum, stable quasi particles can form again. One has 
to describe such a process as with usual pairing but generalised to the case of two 
fermions with unequal massesf® Another possible channel for a condensation is 
its formation out of two deuterons, see Fig. [TUT b). This is similar to pairing of two 
bosons.^' Which of the various processes will in the end win has to be evaluated 
in the future. Another open question is how the quartet order parameter behaves 
as a function of asymmetry and temperature. The dependence on asymmetry can 
partially be deduced from our study of the critical temperature in However, even 
the latter is incomplete, since only valid at the upper end of the density. In the limit 
as the density goes to zero we have to generalise the study of the critical temperature 
in such a way that the one of the ideal Bose gas is obtained. In the pairing case this 
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a) b) 




Fig. 10. (a) Schematic representation of single particle mass operator with a particle condensate 
pairing up a nucleon (n,p) with a trion ( 3 He, t); (b) mass operator for deuteron propagation 
pairing up with another deuteron of opposite momentum. 



has been achieved by Nozieres and Schmitt- Rink.^' It is certainly a challenge to 
incorporate such an approach to quartet condensation. Coulomb repulsion between 
the a-particles may eventually favor an a crystal structure.^ Competition between 
Bose condensed phase and crystalization may be investigated in future. Quartetting 
is not confined to nuclear physics. As mentioned several times, there exists the 
possibility of the formation of a gas of bi-excitons in semiconductors.^ 1 Also, one 
may speculate about the possibility that with cold atoms four different species of 
fermions will be trapped in the future. If they attract one another as the four 
nucleons in nuclear physics, certainly the quartet condensation can also be studied 
in such systems. First theoretical investigations have already appeared. 19 ' The 
theory of quartetting in infinite systems is certainly not complete. More theoretical 
investigations are needed to get a full picture of the process. 
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